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We calculate the shape correlations between trispectra in various equilateral non-Gaussian models, 
including DBI inflation, ghost inflation and Lifshitz scalars, using the full trispectrum as well as the 
reduced trispectum. We find that most theoretical models are distinguishable from the shapes of 
primordial trispectra except for several exceptions where it is difficult to discriminate between the 
models, such as single field DBI inflation and a Lifshitz scalar model. We introduce an estimator 
for the amplitude of the trispectrum, P^l" an d relate it to model parameters in various models. 
Using constraints on gt?^ % from WMAP5, we give constraints on the model parameters. 
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I. INTRODUCTION 



Almost scale-invariant and Gaussian primordial cosmological perturbations predicted by inflation are consistent 
with observational data. While this suggests that inflation did happen in the early universe, there still remain 
many important questions about inflation. One of the most important problems is to distinguish between various 
inflationary models as there are still many models that are consistent with observational data at present. For this 
purpose, it is necessary to have further information about the early universe such as deviations from Gaussianity 
I ■ (non-Gaussianity) of the primordial perturbations and primordial gravitational waves. In this paper, we focus on 
primordial non-Gaussianity. 

Actually, the statistical properties of primordial fluctuations provide crucial information on the physics of the 
very early universe [1-4] (see [5] for a review). In the simplest single field inflation models where the scalar field 
has a canonical kinetic term and quantum fluctuations are generated from the standard Bunch-Davis vacuum, non- 
Gaussianity of the fluctuations is too small to be observed even with future experiments [6-8]. Thus, the detection 
of non-negligible deviations from Gaussianity of primordial fluctuations will have a huge impact on the models of 
the early universe. So far, most of the studies have focused on the leading order non-Gaussianity measured by the 
t^J- \ three-point function of Cosmic Microwave Background (CMB) anisotropics, i.e. the bispectrum [9-11]. Especially, the 
f**«« ■ optimal method of extracting the bispectrum from the CMB data has been sufficiently developed [12-18] (for a more 
general approach, see [19-21]). However, the bispectrum includes only a part of information about non-Gaussianity 
. and many models can predict similar bispectra. 

For example, k-inflation [22, 23] and Dirac-Born-Infcld (DBI) inflation [24] are shown to predict almost equilateral 
type bispectrum (see [25, 26] for reviews). However, future experiments like Planck [27] can also prove the higher 
order statistics such as the trispectrum [28-30] which gives information that cannot be obtained from the bispectrum 
[31-33]. 

Non-Gaussianity of the observed CMB anisotropies comes from not only the primordial origin but also the non- 
linear effects in the CMB at late times such as the coupling between Integrated Sachs- Wolfe (ISW) effect and the 
weak gravitational lensing [34]. However, these non-linear effects can be negligible compared with the primordial 
non-Gaussainity in the models such as DBI inflation and ghost inflation where large primordial non-Gaussianity is 
predicted and it is the dominant contribution to the observed non-Gaussianity. In this paper, we focus on large 
primordial non-Gaussianity. 

While the trispectrum has more information, it requires more work to understand how to distinguish between 
various trispectra predicted in many theoretical models and how to measure the amplitude of the trispectrum because 
it has more parameters than the bispectrum. To estimate an overlap between the shapes of two different trispectra, 
the shape correlator was introduced by Regan et.al [35] based on the reduced trispectrum. Furthermore, based on 
this shape correlator, two of us investigated an estimator to measure the amplitude of the trispectrum in some 

equilateral type non-Gaussian models, like k-inflation, single field DBI inflation and multi-field DBI inflation [36]. 
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Recently, the shape dependence of the trispectra from ghost inflation model [37, 38] and Lifshitz scalar model 
have been calculated in Refs. [39, 40] and Ref. [41]. These models are known to give equilateral type bispectra. For 
Lifshitz scalar model, we will see that there is no natural way to decompose the full trispectrum into the reduced 
trispectra in some cases and the shape correlator based on the reduced trispectrum is not necessarily well-defined. 
Therefore, in this paper, we consider an implementation of the shape correlator of the primordial trispectrum using 
the full trispectrum so that we can calculate the shape correlations and investigate estimators in equilateral type 
non-Gaussian models including ghost inflation and Lifshitz scalar. 

The rest of this paper is organised as follows. In section II, we introduce two shape correlators defined in different 
ways where one is defined based on the reduced trispectrum and the other is based on the full trispectrum. In 
section III, we calculate the shape correlations of the trispectra among various theoretical models; DBI inflation, 
ghost inflation and Lifshitz scalar. In section IV, we express the estimators g^ %l in terms of the model parameters, 
which is useful to constrain the model parameters from future experiments. We also give constraints on model 
parameters using WMAP constraints on g^ tl obtained in [43] . Section V is devoted to the summary and discussions 
of this paper. In appendix sections A, B and C, we review the shape functions of the reduced trispectra in DBI 
inflation model, ghost inflation model and Lifshitz scalar model. 



In this section, we introduce two shape correlators. One is introduced by Regan [35] and we review this method 
in section II A. In this method, the shape correlators are defined based on the reduced trispectrum. Although the 
reduced trispectum includes all information of the original trispectrum, there is no unique way to decompose the 
total trispectrum into the reduced trispectrum. Then there appears an ambiguity in the definition of shape correlator 
based on the reduced trispectrum; a different choice of the reduced trispectrum gives a different shape correlator. 

Before discussing the shape correlators, we review the definition of the trispectrum. The trispectrum 
Tf(ki, k 2 , k 3 , k4) of the curvature perturbation £ is defined as 



where C(ki) is a Fourier component with the momentum ki and the subscript "c" in the left hand side denotes the 
connected component. The trispectrum T^(ki, k 2 , k 3 , k4) generally depends on four three-momenta, namely twelve 
parameters. Assuming isotropy and homogeneity of the universe on large scales, the number of parameters reduces 
to six. 



Here, we review the shape correlator discussed in [35] . First we exploit the symmetry of the trispectrum to define 
the reduced trispectrum as follows [28] . We rewrite the definition of the trispectrum as 



Because of the symmetry of the trispectrum, the reduced trispectrum includes all information of the original trispec- 
trum [28]. However, this decomposition of the trispectrum is not unique because there is an ambiguity in choosing K 
in Eq. (2). In some cases, there is a natural choice of the reduced trispectrum. For example, in the case of trispectrum 
produced by two three-point vertices, it consists of s-, t- and u-channels. Then we can decompose the trispectrum 
into three parts accordingly and define the reduced trispectrum in a natural way. 

The reduced trispectrum 7^(ki, k 2 , k 3 , k 4 ; K) depends on six variables and we can choose them to be 
(k\, fc 2 , /c3, fc 4 , ki2, O4) where O4 represents the deviation of the quadrilateral from planarity which is specified by 
the triangle (k\, k2, &12) (see Fig 1). From the geometric restriction, the range of 64 is constrained as 



II. TWO TYPES OF SHAPE CORRELATORS 



(C(k 1 )C(k 2 )C(k 3 )C(k 4 ))c = (2^) 3 5 3 (k 1 +k 2 +k 3 +k 4 )T c (k 1 ,k 2 ,k3,k 4 ), 



(1) 



A. The shape correlator based on the reduced trispectrum 




(2) 





< COS 64 < 1. 



(3) 



2fci 2 fc 4 



Motivated by the relation between the CMB trispectrum and the trispectrum for £, we define the shape function 



FIG. 1: definition of 4 



for the reduced trispectrum as 1 

Sr{ki,k 2 ,k 3n k 4l ki2,9 i ) = (fcifc 2 fc 3 fc 4 ) 2 fci27c(fci, k 2 , k 3 , fc 4 ; k 12 , 64) ■ (4) 
Regan et al. [35] proposed to define an inner product between two different shape functions Sj- and S77 as 
F(S r ,S T ') 

dkidk 2 dk 3 dkidk 12 / d(cos64)Sr(ki,k 2 ,k 3 , fc 4 , k 12 , 4 )Sy/(fci, fc 2 , fc 3 , fc 4 , k 12 , 6±)w{k\, k 2 , fc 3 , fc 4 , fc 12 ) ,(5) 



where iy is an appropriate weight function. The weight function should be chosen such that S 2 w in k space produces 
the same scaling as the estimator in I space and we adopt the one used in Ref. [35] 2 , 

fci2 

w(k 1 ,k 2 ,k 3 ,k4,k 12 ) = — — — — — , 2 ■ (6) 

(fci +k 2 + k 12 y(k 3 + fc 4 + k 12 y 

The integration range of the momenta fci, k 2 , k 3 , fc 4 , ki 2 in the integral of Eq. (5) is determined by the triangle 
inequality of the momenta. Namely, the conditions to satisfy ki2 = ki + k 2 are 

ki 2 <ki+k 2l ki<ki 2 + k 2 and k 2 < ki 2 + k±. (7) 

Moreover, from the momentum conservation we obtain — ki2 = k 3 + k4, and thus the conditions 

ki 2 <k 3 + ki, k 3 <ki 2 + ki and fc 4 <fci2 + ^3, (8) 

must be also imposed. The integration range of cos0 4 is fixed by inequality (3). Due to the symmetry under 
interchange of fci and k 2 , we can confine the integration range to fci > k 2 . With this choice of weight, the shape 
correlator is defined as 

C(S T , S T < ) = ; F(5r ' 5r,) . (9) 

y/F(ST,ST)F(ST>,S T >) 



1 The factor to relate the reduced trispectrum with the shape function is not unique. For example, in Ref. [35], instead of (A.'i£.'2fc3fc 4 ) 2 fci2, 
another choice (fciA^fcsfc^ 9 / 4 is also discussed. But it is possible to check that the dependence on this factor is not very significant 
when the shape correlation is sufficiently large, like C(S-j-,Sq-i) > 0.7. 

2 The choice of the weight function is not unique, either. For example, in the first version of Ref. [35], instead of , — vifi. it. 11. — ct> 
another choice t — n — , , 1 . ,, — -r-, — -r-. — r is considered. But again, we checked that the dependence on this factor is not significant 

"•12 +^2 + fc 12 ~r^4 +"-12 ) 

when the shape correlation is sufficiently large. 
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B. The shape correlator based on the full trispectrum 

In this subsection, we implement the shape correlator using the full trispectrum Tf . As we mentioned in the previous 
subsection, while the reduced trispectrum has all information of the original trispectrum, the decomposition of the 
full trispectrum into the reduced trispectra is not unique. If all momenta appear symmetrically in the definition of the 
correlator, the obtained correlation is independent of the choice of the reduced trispectrum. However, the definition 
given by Eq. (5) apparently breaks the symmetry because of the integration variables k\2 and 64. Therefore, in order 
to define the correlator uniquely, we must use the full trispectrum in the definition of the shape correlator. 

The shape correlator based on the full trispectrum is defined in a similar way with the one based on the reduced 
trispectrum. Namely, we begin with defining the shape function as 

S T (ki,k2,k 3 , k 4 , ki2,0 4 ) = (k 1 k 2 k 3 k4) 2 k 12 T c (k 1 , fc 2 , fc 3 , fc 4 ; fc 12 , 64) • (10) 

The inner product is defined as 

F(St, St<) 

dkidk 2 dk 3 dk4dki2 J d(cos 04)<Sr(fci, k 2 ,k 3 ,k4, k-42, 04,)ST'(ki, fc 2 , fc 3 , k 4 , k-42, 9i)w{ki, fc 2 , fc 3 , fc 4 , fci 2 ) .(11) 

where we use the same weight function w as that in the case of reduced trispectrum. The domain of integration in 
Eq. (11) is basically the same as that for the reduced trispectrum, i.e. inequalities (3), (7) and (8). However, since 
the symmetry under interchange of k\ and fc 2 still remains, the additional condition ki > fc 2 does not change the final 
value of the shape correlator. Therefore, we adopt inequalities (3), (7), (8) and k\ > fc 2 as the domain of integration. 
The shape correlator is defined in the same way 

C{St,St>) = , F{St ' St,) . (12) 

\J F(St, St)F(St> , St 1 ) 



III. SHAPE CORRELATIONS 



In this section, we explicitly calculate the shape correlations among trispectra predicted by various theoretical 
models. The main purpose of investing the trispectrum is discriminating between models which are hard to be distin- 
guished by the bispectrum. Thus we concentrate on the shape correlations among the models where the bispectrum 
is dominated by the equilateral type one; single field and multi-field DBI inflation models, ghost inflation model and 
Lifshitz scalar model. 

The explicit forms of the shape functions are calculated in Appendixes. S® BI ^ a \ S BBI ^ and S^ ost are defined as 
the shape functions in single field DBI Inflation, multi-field DBI inflation and ghost inflation, respectively. Explicit 
forms of them are written in Eqs.(A32), (A35) and (B3). In Lifshitz scalar field, there are a few contributions of 
Trispectrum and it depends on the parameters in the theory. We define S^ se ' U \ S^ se ' 12 \ c^_( se ' 13 ) ^ S^ cl,1 \ st^ cl ' 2 \ 
and as the shape functions of respective contributions. We show explicit forms of them in Eqs.(C2) and (C3). 

Furthermore, the shape given by (see Eq. (A10)) can be written in a separable form [42] and it provides a fast 
estimator for the equilateral type trispectrum [35, 36]. This means that if the shape function for the trispectrum is 
highly correlated with S7J , we can constrain the model parameters from the amplitude of the trispectrum in a very 
simple and fast way as we will discuss in Sec. IV. Therefore, we also calculate the shape correlations of the trispectra 
predicted by theories mentioned above with 

Another motivation for calculating the shape correlations is to check the difference between the two shape correlators 
that we introduced in the previous section. As we mentioned in the previous section, the shape correlation based on 
the reduced trispectrum depends on the way to decompose the full trispectrum into the reduced trispectra. Thus 
strictly speaking, we should use the full tripsectrum when calculating the shape correlation. However, the calculation 
is easier for the shape correlation using the reduced trispectrum. If the difference is insignificant, we may still use the 
reduced trispectrum to calculate the shape correlator. 



A. Shape correlations based on the reduced trispectrum 



Here, we calculate the shape correlations among the theoretical models based on the reduced trispectrum. We must 
specify how we decompose the full trispectrum into reduced trispectra. One natural way is to decompose it so that 



FIG. 2: scalar exchange diagram with ki2 inner propagator 



one of the reduced trispectrum depends only on five parameters fci, fc 2 , k%, k^ and fci 2 . Except for the case of the 
trispectrum obtained by the contact interaction in Lifshitz scalar model for i = 3 in Eq. (C3), it is possible to find 
such a decomposition. In fact, in the case of the scalar exchange trispectrum, this is always the case. The contribution 
of the trispectrum from the scalar exchange can be decomposed into s-, t- and u-channels and one of the channel 
depends only on five parameters hi, fc 2 , k 3 , fc 4 and fc 12 . This can be understood as follows. The diagram of the scalar 
exchange trispectrum has one inner propagator. In the diagram where the momentum of the inner propagator is 
fci2, the propagators with momenta fci and fc 2 meet directly at one vertex. On the other hand, the propagators with 
momenta ks and k^ meet at the other vertex (see Fig. 2). Then, the trispectrum from Fig. 2 is written in terms of 
fci, &2, fe, ki, k\2, ki • k2 and k3 • k4. The inner products ki • k2 and k3 • ki can be expressed as (fc 2 2 — k\ — fc 2 )/2 
and (fc 2 2 — fcf — k%)/2, respectively. Therefore, the trispectrum from Fig. 2 depends only on fci, fc 2 , fc3, k± and fc i2 . 
Moreover, in most models even the contact interaction trispectrum can be also decomposed in such a way that the 
reduced trispectrum only depends on fci, fc 2 , ^3, ki and fci 2 . In the contact interaction, there are no inner propagators. 
Thus the dependence on fci 2 , fc i3 and fci4 stem only from the derivative coupling such as (Vi</>) 2 </> 2 . Only when the 
product of fci 2 and fci3 appear in a contact interaction trispectrum, the trispectrum can not be decomposed in the way 
so that reduced trispectrum depends only on fci, fc 2 , fc3, k^ and fci 2 . This actually happens for the contact interaction 
in Lifshitz scalar field for i = 3 in Eq. (C3). 

In appendix A, B, C, we summarise the explicit expressions of the shape functions defined by the reduced trispectra 
in single field and multi-field DBI inflation models, ghost inflation model and Lifshitz scalar model, respectively. As is 
mentioned above, except for the one corresponding to the contact interaction in Lifshitz scalar field for i = 3 in (C3), 
they depend only on five parameters fci , fc 2 , k 3 , k^ and fci 2 . For the trispectrum coming from the contact interaction 
in Lifshitz scalar field for i = 3 in Eq. (C3), we need to define the reduced trispectrum so that it is symmetric under 
the transpose among fci, fc 2 , k a and k 4 . 

The shape correlations among the models are summarised in Table I. From this table, S® BI ^ and 5*^ SL ' 22 ) t urn 
out to be highly correlated with S7- 1 . 



Now, we move to the shape correlations based on the full trispectrum. With Eqs. (1) and (2), the full trispectrum is 
obtained from the reduced trispectrum summarised in the appendix sections, by adding their permutations (fc 2 -f* ks) 
and (fc 2 -f-s- k 4 )- Then, the shape functions based on the full trispectrum can be written by fci, fc 2 , k^, k 4 , fci 2 , fci3 and 
fci4. We substitute the shape functions of the full trispectra into the definition of the shape correlator (12). Thus, we 
need to express fci3 and fci4 in terms of fci, fc 2 , fc3, fc4, fci 2 and 64. We find that these variables are expressed as 



B. Shape correlations based on the full trispectrum 




1 



12 



^4fc 2 fc 2 2 - (fc 2 + fc 2 2 - k\f y / 4fc|fc 2 2 cos 2 9 A - (fc 2 + fc 2 2 - fc 2 ) 2 ) 





1 



12 



^k\k\ 2 - (fc 2 + fc 2 2 - k 2 2 f^k\k\ 2 cos 2 6> 4 - (fci + fc 2 2 - fc|) 2 ) 




(13) 



where ± is determined by the configurations of the momenta; if the angle between the planes which are specified by 
(fci, fc 2 , fci 2 )-triangle and (^3, k^, fci 2 )-triangle, a, is smaller than the right angle (see Fig. 3), the corresponding sign is 
the minus (— ). Otherwise, we should choose the plus (+). 
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TABLE I: The shape correlations among the models based on the reduced trispectra 







Q DBI(s) 


qghost 


n/i(se, 11) 


Q h(se, 12) 


Q h(se,22) 




Mci,2) 


Qh(ci,3) 




0.87 


0.33 


0.24 


0.24 


-0.62 


0.95 


-0.35 


0.01 


-0.53 


gDBI(cr) 




0.19 


0.41 


0.42 


-0.78 


0.96 


-0.60 


0.05 


-0.67 


gDBI(s) 






0.59 


0.13 


-0.25 


0.28 


-0.21 


0.20 


-0.78 


Qtghost 








0.23 


-0.31 


0.30 


- 0.86 


0.32 


0.36 


Q,h(se,ll) 










-0.75 


0.33 


-0.37 


0.05 


-0.25 


Q h(se,12) 












-0.73 


0.50 


-0.07 


0.52 


rih(se,22) 














-0.44 


0.02 


-0.63 


















-0.25 


0.19 


gh(ci,2) 


















0.08 



TABLE II: The shape correlations among the models based on the full trispectra 





gDBI(a) 


qOBI(s) 


qghost 


gh(se, 11) 


gh(se, 12) 


gfc(se,22) 


gh(ci,l) 


gh(ci,2) 




QCl 


0.96 


0.42 


0.41 


0.23 


-0.70 


0.99 


-0.54 


-0.01 


-0.38 


gDBI(,T) 




0.40 


0.51 


0.36 


-0.80 


0.98 


-0.71 


-0.01 


-0.43 


Q DBI(s) 






0.76 


0.25 


-0.45 


0.45 


-0.39 


0.16 


0.40 


Qtghost 








0.24 


-0.42 


0.48 


-0.83 


0.16 


0.10 


^h(se, 11) 










- 0.70 


0.32 


-0.36 


0.02 


-0.06 


gh(se,12) 












-0.77 


0.60 


-0.01 


0.22 


gh(se,22) 














-0.62 


-0.01 


-0.38 


qh(ci,l) 
















-0.07 


0.28 


sr> 2) 


















0.07 



The shape correlations among the models are summarised in Table II. In most of the models, the shape correlations 
based on the reduced trispectra are close to those based on the full trispectra. Especially, if the correlation based 
on the reduced trispectra is high, there is a little difference between the two shape correlations. Thus as long as the 
shape correlation is high, it is possible to use the reduced trispectrum to calculate the shape correlation. On the other 
hand, a care must be taken if the correlation is low. 



IV. THEORETICAL PREDICTIONS AND OBSERVATIONAL CONSTRAINTS 

In this section, in order to constrain the model parameters from the amplitude of the trispectrum, we will calculate 
the amplitude of the estimators <7nl" based on the shape (see Eq. (A10)), which can be written in a separable 
form [42]. We also discuss the amplitude of the trispectrum in the regular tetrahedron limit, iNL and the maximally 
symmetric configuration, tnl- 



A. g e N T l 

Here, following Ref. [36], we discuss the theoretical predictions for the amplitude of the trispectrum ■ We first 



define the equilateral shape using S? 1 



T ' 



gepuil = 6^2^)3^ (14) 



7 




FIG. 3: definition of a 



where is defined as (^(ki^k^)) = 167r 5 <5 3 (ki +k2)fc 1 3 V(. The parameter is defined as the normalised shape 

correlation between Sj- and S^ ml ; 

Tp (a Qequil\ 
equil _ f l J T,Jr > 
"NL j-^/^equil riequil \ ' v ' 

i<(b T ,b r ) 

We explicitly calculate t^L %l m single field DBI inflation model, multi-field DBI inflation model, ghost inflation 
model and Lifshitz scalar model. The results are summarised in Table III. Since we have two kinds of the shape 
correlators (one is based on the reduced trispectrum and the other is based on the full trispectrum) , we show the two 
results in Table III. Here, we use the fact that the power specrta of the primordial perturbations in single field and 
multi-field DBI inflation models, ghost inflation model and Lifshitz scalar model are given by 



1 H' 
2^2 2^2' 

1 H a T-jis 

1 f H 



Pf = ^^- 1171 



V " h ° St ^(r(l/4)) 2 a3/4 [m) ' ,LS) 



' c 27T 2 2/i 2 - 



Moreover, in Lifshitz scalar model, we individually show the contributions from each shape function. 

Notice that although we obtain g^ 11 for all theoretical models mentioned above, this <7nl" can ^ e meaningful 
only in the case where the shape correlation with the equilateral shape is sufficiently high, like > 0.7. Otherwise, this 
template based on the equilateral shape S^ ul1 does not fit the shape of the trispectrum well and the noise dominates 
over the signal 3 . 

From this table, we find that there are little differences between g^f based on the reduced trispectrum and that 
based on the full trispectrum in most of shape functions, especially for S^- BI , ghi se > 12 ) an( j gft(se,22) w j lere 
shape correlation with S^ utl is high enough and the use of <?nl" to measure the amplitude is justified. 



3 We thank D. Regan for pointing out this point. 
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TABLE III: The value of based on reduced trispectrum and full trispectrum. 



DBI(tr) | gDBI(s) 



-tghost 



ih(se,U) 



-rh(se,12) 



reduced 



1.4 x 104 



2.3 



r 4 T 2 
c s 1 -R,S 



4.9 x 10 



5 _5_ 



5 I 4.8x10-5 



1.6 x 10 7 6x 



-1.8 x 10 8 6i& 2 



full 



1.3 x 104- 



2-7^ 



6.0 x 10 



5 7 



,1/2 



875 I 4.8X10" 



1.1 x W 7 bl 



,1/2 



-1.6 x 10»& 1 fe 2 



,1/2 



^h(se,22) 



~ih(ci, 1) 



h(ci,2) 



s: 



h(ci,3) 



reduced 



4.5 x 10 8 b| , 4 . 8X10 _ 



-9.1 x 10'ti 



v 1/2 



7.5 x 10 7 £ 2 



-3.7 x 10 8 t 3 



4.8X10" 5 



full 



4.3 x 10 8 b| 



-9.3 x 10 7 ti 



-3.7 x 10 7 i 2 



,1/2 



-6.2 x 10 7 t 3 ' — - 



,1/2 



4.8X10" 5 



B. i NL and tatl 

In order to estimate the amplitude of the trispectrum with different shapes, Chen ct al. [42] define the amplitude 
of the trispectrum £nl using a particular configuration as follows; 

(C 4 )-^(2^) 9 7' c 3 ( 5 3 (^k i )lt NL , (20) 

i 

where the limit stands for the regular tetrahedron limit (fci = ^2 = ^3 = ^4 = ^12 = &13 = k). Theoretical predictions 
for £nl are summarised in Tables IV. 

It is clear that £nl gives only a rough estimation of the amplitude of the trispectrum as it is obtained by using only 
one specific configuration. In fact, for S^ cl ' 3 \ £nl is because S^ 1 '^ happens to vanish for the configuration given 
by (ki = k,2 = k$ = ki = k\2 = &13 = k). On the other hand, if we consider all the configurations, there is a small 
but non- negligible over-lap between this shape and S^. This clearly demonstrates that it is necessary to use all the 
configurations to estimate the amplitude of the trispectrum. 

We also show the amplitude of the trispectrum in another definition. The configuration {k\ = &2 = k$ = ki = k 
and fci2 = fci3 = fci4 = 2fc/v / 3) is maximally symmetric while the configuration of tetrahedron is not invariant 
under interchange of ki and k 2 . Therefore, the maximally symmetric configuration is probably useful to analyze the 
trispectrum. Here, we define the amplitude of the trispectrum tnl as in [56]; 4 

TNL = ir/wk?' (21) 

where is the trispectrum of the maximally symmetric configuration. Theoretial predictions for t^l are also 
summarised in Tables IV. 



C. Constraints on model parameters 

Observational constraints on the equilateral non-Gaussianity were obtained from WMAP5 data by applying the 
estimator with the shape given by [43]. The constraints on the amplitude of the trispectrum in the regular 
tetrahedron limit is obtained as 

t equil = ( _ 3 n ± ? 5) x 1Q 6 ( 68%CL ). (22) 

This constraint can be converted to that on as follows. We define the following shape function 

64 

S r = 9 c Z%(^ 2 Vt?S?, (23) 



4 While in [56] the amplitude of the trispectrum tnl is defined in any configuration, we rewrite the simplified form which can be applied 
only in the maximally symmetric case. 
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TABLE IV: The value of £nl and t nl . 



qDBI(ct) 



~,DBI(s) 



Qghos 



R7T 



Mse,ll) 



-,h(se,12) 



5.4 x io -1 ^ 



1.4 x lO- 1 -^ 



-1.2 x 10 5 -i» 



4.8X10- 5 



5.1 x 10 6 6f 



-2.3 x 10 7 &ifc 2 



TNL 



5.6 x lcr 1 ^ 



1.2 x 10 _1 



-1.0 x 10 



,4 7 



5 I 4.8x10- 



1.2 x 10 7 b? 



4.8XK)- 5 



-9.9 x 10 7 6i&2 



s: 



h(se,22) 



-»/i(ci,l) 



-,h(ci,2) 



Qh(ci,3) 



TNL 



2.8 X 10 7 &2 4 



2.7 x 10% 



4.8X10- 5 



6.9 x 10 J t 2 



4.8X10- 5 



2-1 x 10 8 bi ( 



4.8X10- 5 



3.1 x 10% ( 21B ^ r 



,1/2 



3.1 x 10% ' — ^ 



,1/2 



4.8X10- 5 



TABLE V: Constraints on model parameters from g^L**' m various models. 



reduced 



full 



~,DBI(a) 



1.8 X 10 _J < C s 



1.7 x 1Q~ Z < c s 



~,DBI(s) 



1.1 x lO" 2 < CsVTtzs 



1.2 x KT 2 < CsVTns 



^75" 



Qghost 



,!/2 



6.9 x 10 2 < ^ TJ ^ FW < 2.9 x 10 



-5.7 x 10' < ^ I 



< 2.3 x 10" 



-,h(se,ll) 



(0 <) bi 



r c /2 

1 \ 4.8X1C)- 5 



< 



(0 <) b 



2 1 

1 I 4.8X 10~ 5 



< 1.3 x 10 



-ih(se,12) 



-7.8 X lO" 1 < &1&2 



pV2 



3.8 x 10- 1 < 6162 



< 2.1 



-,h(se,22) 



(0 <) 62 I T7-7T7- 



< 3.1 x 10" 1 



(° <) & 2 ( 4"5xW ) < 3 3 x 10_1 



-~,h(ci,l) 



v 1/2 

-1.5 < ti I , . < _ t ) <3.7 



-1.5 < ti 



,V2 



< 3.7 



-,h(ci,2) 



-4.5 < * 2 



,V2 



< 1.9 



-3.8 < t 2 



-3.8 x IP' 1 <t 3 ( j£o. t ) < 9.2 x HT 1 



< 9.2 



-,h(ci,3) 



-2.3 < t 3 ( ] < 5.5 



which by definition gives wnen we apply Eq. (15). Now using the relation between the shape function and the 

full trispectrum and evaluating the full trispectrum in the regular tetrahedron limit, we can calculate t^ lt for this 
trispectrum as 

,equil 1 equil (OA\ 

nl — 32 5nl ' ^ ' 

Then we obtain the constraints on <7nl" as 

g equil = ( _ 9 95 ± 24) x 1Q 7 ( 68%CL )_ (25) 

From the theoretical predictions for <?nl" in various models obtained by using the full trispectrum as well as the 
reduced trispectrum, we can derive the constraints on model parameters as in Table V. 



V. SUMMARY AND DISCUSSIONS 



There are many interesting early universe models motivated by string theory and the effective field theory that 
predict large primordial non-Gaussianity with the equilateral type bispectrum. Given that future experiments like 
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Planck can prove even higher-order statistics, it is important to investigate whether the shape of the trispectrum can 
distinguish such equilateral type non-Gaussian models. For this purpose, the shape correlator of the trispectrum is 
particularly useful. 

So far, the shape correlator constructed from the reduced trispectrum has been adopted [35]. While the reduced 
trispectrum has all information of the full trispectrum, the shape correlator based on the reduced trispectrum depends 
on the way in which the full trispectrum is decomposed into reduced trispectra because the form of the integration 
variables in the shape correlator breaks the transposition invariance of momenta ki, k 2 , k 3 and k 4 . In most equilateral 
type non-Gaussian models such as DBI inflation model, there is a natural way to decompose the full trispectrum into 
the reduced trispectra so that one of the reduced trispectra depends only on five parameters fci, k%, k^, k^ and ki2- 
However, for some classes of equilateral type non-Gaussian model like Lifshitz scalar model, this decomposition is not 
possible. Therefore, in this paper, we studied the shape correlator of the primordial trispectrum based on the full 
trispectrum. 

In order to check the difference between the two shape correlators, we calculated the shape correlations among 
trispectra in various equilateral non-Gaussian models; DBI inflation, ghost inflation and Lifshitz scalar models. We 
found that both shape correlators give similar results as long as the shape correlations are high. 

From the shape correlations, it is possible to judge whether we can distinguish between various equilateral non- 
Gaussian models using the tripsectrum. For example, we showed that it is difficult to distinguish between the single 
field DBI inflation model and the 62&2-dominated Lifshitz scalar model by the shape of the trispectrum. Since both 
single field DBI inflation model [44] and &2^2-dominated Lifshitz scalar model [41] predict almost the same equilateral 
type bispectrum, we can not distinguish between these two models by the primordial non-Gaussianity up to this order. 
In order to distinguish between these models, we need other information such as the primordial gravitational wave. 
A similar conclusion holds for the comparison between the ghost inflation model and ^-dominated Lifshitz scalar 
model. Despite these exceptions, our result suggests that we can distinguish between many equilateral non-Gaussian 
models, which predict almost the same equilateral type bispectrum, from the shape of the trispectrum. 

On the other hand, in order to measure the amplitude of the trispectrum and constrain parameters in the theoretical 
models, it is necessary to develop an estimator for the trispectrum. Since the form of the trispectrum is too complicated 
for this class of models, it is generally impossible to construct an optimal estimator (see however [35] for the model 
independent approach). Using the fact that the shape given by Eq. (A10) can be written as a separable form, it is 
possible to construct a fast optimal estimator <?nl" [36j 43]. We expressed the estimator in terms of models parameters 
in various models. The constraint on <7nl" was obtained from WMAP5 in [43]. From this constraint, we obtained 
constraints on model parameters. We emphasized that the amplitude of the trispectrum for a particular configuration 
such as the regular tetrahedron limit cannot be reliably used to characterise the amplitude of the trispectrum and we 
need to use all configurations to define the amplitude <7nl" ■ 

Finally, it is known that that inflation models based on Galileon and its generalisations, which are shown to be 
the most general single field inflation model with second-order field equations [45, 46] give also the equilateral type 
bispectrum [47-54] (see also for the discussion about the shape dependence of the bispectrum in this type of inflation 
model [55]). It would be interesting to study the amplitude of the trispectrum j^" 1 in these models. 
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Here, based on our previous work [56] , we summarise the shape functions for the reduced trispectra in general single 
field k-inflation models described by the following action: 
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Appendix A: Shape functions in general single field k-inflation 




(Al) 
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where (j> is the inflaton field, R is the Ricci scalar and X = — (l/2)g^d^(j)d v (l>, where g^ u is the metric tensor. 

In this class of models, the third and the fourth order interaction Hamiltonian of the field perturbation S(f> in the 
flat gauge at leading order in the slow-roll expansion are given by 



I 



d 3 x 
d 3 x 



Aa5(f)' 3 + Bade/)' (ddc/)) 2 



(A2) 
(A3) 



where prime denotes derivative with respect to conformal time rj and coefficients A, B, /3i, /? 2 and [}$ are given by 



<2X 



P,xx + T^XP.XXX 



B 



f2X 



,xx- 



(A4) 



= P 



Pi 

P2 = 
P 3 = 



XX 



-\p,xx 



2XP XXX 
1 



4 s 



1 



X 3 C 2 S 2 

~5~ Lp ,xxx - «A^P 4 x, 



+-XciP X xx, 



c 2 

^P 



,xx- 



(A5) 



Here, Pj denotes the derivative of P with respect to X, P t xx denotes the second derivative of P with respect to 
X, and so on. c s is the sound speed of the perturbation of the scalar field <j) which is defined as 



:x 



Px + 2XP, 



The shape function Sj, is composed of two parts 



ok _ nk(cont) 



,xx 



k(scalar) 
T 



(A6) 



(A7) 



where s^ cont ^> denotes the contribution from the contact interaction and g^ bcaia ' > denotes that from the scalar 
exchange interaction, respectively. 

The shape function for the reduced trispectrum arising from the contact interaction g^ cont ^ depend on five param- 
eters fei, fc 2 , fc 3 , hi and k\2- It is given by 



-ik(scalar) 



s Hcon t) = (_ 24j8lC 3 S c, _ p 2CsS c, _ 2/33Cs -1^3) 1_ N S . 



Here , and arc the following shape functions: 
k 12 Uf =1 ki 



(El^)' 



+ 3 perms. 



(A8) 



(A9) 



fci 2 fc 2 fc|(k 3 ■ k 4 ) 

L(Eti^) 3 nf =1 fci 

fc^fcffcKk! • k 2 ) 



+ 



Sf 3 



/ci 2 (k! • k 2 )(k 3 • k 4 ) 




+ 3 perms. . 



Ei<j 



12- 



n 



(A10) 



3 perms. , (All) 



where "3 perms." denotes the permutations (ki -H- k 2 ), (k^ -H- £4) and (fci -H- £;2,fc3 «-> £4). In Eq. (A8), _ff is the 
Hubble parameter at inflation era and N = H/ ^JIP xc-s ■ 
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Similarly, S. 



k(scalar) 
T 



is given by 



^(scalar) 



{A 2 c 4 s S^ + ABc 2 s S¥ + B 2 S s j) 



~ 2 H 2 N 10 
8X 2 



(A12) 



Here S!^ , and S!j? are the following shape functions: 



Sf 1 



-9/ci 2 (fclfc2&3&4) 



1/2 



Ji(fci, fc 2 , -fci2, fc 3 , fc 4 , fci 2 ) - J 7 i(-k 1 ,-k 2 , -fci2, fc 3 , fc 4 , fci 2 ) 



+7"l(fc 3 , fc 4 , -fci2, fcl, fc 2 , fci 2 ) - Ji(-fc 3 , -fc 4 , -^12, fc3, &4, fcl2) 



+ 3 perms. . 



(A13) 



Sf 2 
Sf 2 " 



Sf 26 



5f 2d 



QS2o I 0^26 I QS2c I C s 2d 

-fcl2(fclfc 2 fc 3 fc 4 ) 1 / 2 (k 1 .k 2 )(k 3 -k 4 ) 



(A14) 



J r 2 (-k 12 ,k 1 ,k 2 , fci2, fc 3 , fc 4 ) - T 2 {-k X2 , -fci, -fc 2 , fci2, fc 3 , fc 4 ) 



+^2(-fcl2, fc3, fc4, fcl2, fcl, fe) - ^2(-fcl2, ~fc 3 , "^4, &12, fcl, fe) 



+ 3 perms. 



(A15) 



= -2fc 12 (fcifc 2 fc 3 fc4) 1/2 (ki • k 2 )(k 12 • k 4 ) 



J" 2 (-fcl2, fcl, ^2, fc 3 , fc4, fcl2) - ^2(-fcl2, -fcl, -fc2, fc3, fc4, fcl2) 



+J" 2 (fc 3 , fc4, -fcl2, fcl2, fcl, fc2) - ^2(-fc3, -fc4, -&12, &12, fcl, fc 2 ) 



+ 3 perms. . 



(A16) 



= 2fci 2 (fc 1 fc 2 fc 3 fc 4 ) 1 /2( kl2 . k2 )( k3 . k4 ) 



J" 2 (fcl, fc2, -fcl2, fcl2, fc3, fc 4 ) ~ ^2 (-^1, -^2, -fcl2, fcl2, fc 3 , fcl) 



+J" 2 (-fcl2, fc 3 , fc4, fcl, fc2, fcl2) - ^2(-fcl2, -fc 3 , -fc4, fcl, fc2, fcl2) 



+ 3 perms. . 



(A17) 



4fci2(fcifc 2 fc 3 fc4) 1/2 (ki2 • k 2 )(k 12 • k 4 ) 



J" 2 (fcl, fc2, -fcl2, fc 3 , fc4, fcl2) - ^2(-fcl, -fc2, -fcl2, fc 3 , fc4, fcl2) 



+J" 2 (fc 3 , fc 4 , ~fcl2, h,k 2 ,k 12 ) - J" 2 (-fc 3 , -fc4, -fcl2, fcl, fc2, fcl2) 



+ 3 perms. 



(A18) 



S S T 3 



QS3a I QS3b _|_ OS3c _|_ C«3d 

3fci2(fcifc 2 fc 3 fc 4 ) 1/2 (k 3 -k 4 ) 



Tz{k\,k 2 , -fci2, fci2, fc 3 , fc 4 ) - J" 3 (-fci, -fc 2 , -fcl2, fcl2, fc 3 , fc4) 



+J 4 (-fcl2, fc 3 , fc 4 , fcl, fc 2 , fcl2) - ^4(-fcl2, ~fc 3 , ~fc4, fcl, fc2, fcl2) 



+ 3 perms. . 



6fci 2 (fcifc 2 fc 3 fc4) 1/2 (k 12 ' K 4 j 



•F 3 (fci, fc2, -fci2, fc 3 , fc 4 , fci 2 ) - J" 3 (-fci, -fc 2 , -fci2, fc 3 , fc 4 , fcl2) 



+J" 4 (fc 3 , fc 4 , -fci2, fcl, fc 2 , fcl 2 ) - J" 4 (-fc 3 , -fc 4 , -fcl 2 , fcl, fc 2 , fcl 2 ) 



+ 3 perms. . 



3fc 12 (fcifc 2 fc 3 fc 4 ) 1/2 (ki-k 2 ) 



^4(-fcl2, fcl, fc2, fc 3 , fc 4 , fcl2) - ^4(-fcl2, -fcl, -fc2, fc 3 , fc 4 , fcl2) 



+J" 3 (fc 3 , fc 4 , -fci2, k 12l ki,k 2 ) - J" 3 (-fc 3 , -fc 4 , -fcl2, fcl2, fcl, fc2) 



+ 3 perms. . 



= -6fci 2 (fcifc 2 fc 3 fc 4 ) 1/2 (k 12 -k 2 ) 



J" 4 (fcl, fc 2 , -fcl2, fc 3 , fc 4 , fcl2) - J" 4 (-fcl, -fc2, -fcl2, fc 3 , fc4, fcl2) 



+J" 3 (fc 3 , fc 4 , -fci2, fcl, fc 2 , fcl 2 ) - J" 3 (-fc 3 , -fc 4 , -fcl 2 , fcl, fc 2 , fcl2) 



+ 3 perms. , 



(A19) 



(A20) 



(A21) 



(A22) 



(A23) 



where again "3 perms." denotes the permutations (fci o fc 3 ), (fc 3 <H> fc 4 ) and (fci -R- fc 2 ,fc 3 <H- fc 4 ). Here we have 
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defined four T% functions (with i = 1, . . . , 4) as follows; 
Ti(ki,k 2 ,k 3l k il k 5 ,k 6 ) = -A\kik 2 k 3 k4k 5 k e 
^2(ki,k 2 ,k 3 ,k 4 ,k 5 ,k 6 ) -- 



1 L„i / 



(A24) 



\k 2 k 3 k 5 k 6 \2 AC 



^ k 5 + kg 2 _fc 5 fc 6 



.4 



,4 2 



+ \ [k 2 + fc 3 + fc 5 + fc 6 + I ((fc 2 + fc 3 ) (fc 5 + fc 6 ) + 2k 5 k 6 ) + 2 hh ^ 2 +ka) 

2 , 1 
+ -^(k 5 h + (k 2 + k 3 ) (k 5 + k 6 ) + k 2 k 3 + — (k 2 k 3 (k 5 + k 6 ) + 2k 5 k 6 (k 2 + k 3 )) 



+2 



C 2 

k 2 k 3 k 5 k 6 , 6 



^ 3 (fci,fc 2 ,fc 3 ,fc 4 ,fc 5 ,fc 6 ) 
J r 4{ki,k 2 ,k 3 ,k4,k 5 ,k 6 ) 



+24 
2 



A 2 ' C 3 
k 2 k 3 k 5 k 6 



T^t k 2 k 3 (k 5 + fc 6 ) + k 5 k 6 (k 2 + k 3 ) + 2 



k 2 k 3 k 5 k 6 \ 
A ) 



C 4 

fclfc 2 fc 3 fc 4 |5 1 



|fc 5 fc 6 



AC a 



|A;ifc 4 fc 5 fc 6 |2 1 



1^2 ^3 



^ 3 C 



" fcs + fcg fc 5 fc 6 3 / fc 5 fc 6 \ hk 6 

1 + — — +2-^ + - ^5 + ^ + 2^+12-^- 

1 ^ .4* fc 2 + fc 3 o ^(fc 2 + fc 3 ) + fc 2 fc 3 

+ c + c 2 + c + c 2 

+3^ (^(fc 2 + fca) + 2fc 2 fc 3 ) + 12fefc 3 ^ 



(A25) 
,(A26) 

(A27) 



where .A is defined by the sum of the last three arguments of the Ti functions as A = fc 4 + k§ + k$ and C is defined 
by the sum of all the arguments as C = k\ + k 2 + k 3 + fc 4 + k$ + k§ . 

To simplify the calculation, it is useful to notice that the following properties 



F(5fS S T ) 
F(S t , S t ) 



^a t F(S c T \ S T ), 

i 

Y,a iaj F{S T , Sir).. 



(A28) 
(A29) 



hold for the shape function given by 



s r = J2 a * s: 



r ' 



(A30) 



where i = ci, c 2 , c 3 , Si, s 2 , s 3 and a^s are corresponding coefficients. 

Especially, in the case of single field DBI inflation, the coefficients in the Hamiltonians (A2) and (A3) are given by 



A = — 



1 



ft 



2^1 
1 



B = 



1 



3 ' 



20c; 



2c^ 2 ' ' " 4cf0 2 ' ' " 8c s ^> 2 
and then the shape function based on the reduced trispectrum becomes 



H 



12 



64 



64 



64 



r 



(A31) 



(A32) 



In multi-field DBI inflation model [57], in addition to the shape functions S^, we find it 

convenient to define the following shape functions S^ and S^ given by 



QS2 QS2a QS2b QS2c I QS2d 

O-y- I I I I ' 



C s 3a Q s 36 I QS 3c 

&T T T 



Sf< d 



(A33) 
(A34) 
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With the above functions, the shape function based on the reduced trispectrum of the multi-field DBI inflation model 
can be expressed as 



a DBI{s) ti 



-Is. 



1 



1 



1 



r 



— 5f + — s ¥ + ~ s - 



576 



64 



r 



192 



(A35) 



where T-jis is the transfer coefficient that relate the amplitude of original entropy perturbations to the final curvature 
perturbation 5 . 



Appendix B: Shape function in ghost inflation 



In this appendix, based on our paper [39], we review the trispectrum in ghost inflation and summarise its shape 
functions. Ghost inflation is an inflation model where the inflation is driven by a scalar field <j> in the ghost con- 
densation model. The ghost condensation is the simplest Higgs phase for gravity in infrared and in this model the 
four dimensional diffcomorphism is spontaneously broken by the timelike vacuum expectation value of the derivative 
of the scalar field (f>. Thus, the action for the perturbative field n of <fi does not invariant under four dimensional 
diffcomorphism and relevant terms up to fourth order are written generally as 



/ 



dtdx 3 a 3 



2M 2 



d t TT 5— 



2M 2 



7 (Vtt) 4 
8M 4 o 4 



(Bl) 



where a, (5 and 7 are dimcnsionless constants of order unity and M is a constant with a dimension of mass. Then, 
interaction Hamiltonian is written as 



dtd 3 xa 3 



2M 2 



8M 4 



(B2) 



where 7 = 7 + 2/3 2 . 

Trispectrum from the tree level contribution can be decomposed into two parts. One is obtained by the fourth 
order interaction Hamiltonian and proportional to 7 which is called the contact interaction contribution. The other 
is obtained by the product of the third order interaction Hamiltonian and independent of 7 which is called the scalar 
exchange contribution. 

Here, for simplicity we work only on the contact interaction contribution of the trispectrum as only this contribution 
has new information related with the four-point vertex. From Eq. (B2), if the condition 7 3> f3 2 is satisfied, this 
treatment can be justified. 

According to [39] , the shape function of the contact interaction contribution from the ghost inflation based on the 
reduced trispectrum can be written as 



-tghost 

'r 



7 



H 



2 3 a 3 / 2 \M J V r (V4) 

pO 



7T 



1/2 

(k 1 k 2 k 3 k 4 ) 1 fci 2 (ki • k 2 )(k 3 • k 4 ) 



((-v') 3/2 Hi%( q W 2 )) ((-v') 3/2 H^( q W 2 ))}, (B3) 

where % = y/aH kf/(2M). Here, we choose the reduced trispectrum such that depends on the 5 parameters hi, fc 2 , 
^3, ki and ku as is the case in DBI inflation (see Appendix A). 



5 Strictly speaking, there is another contribution to the trispectrum in multi-field DBI inflation as pointed in Ref. [58]. Since the degree 
of the component studied in [58] depends on background dynamics strongly, we concentrate on the contribution coming from the 
intrinsically quantum four-point function, for simplicity. 
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Appendix C: Shape function in Lifshitz scalar 

In this appendix, based on our previous paper [41], we derive the shape function for the trispectrum in Lifshitz 
scalar. Generally, some shapes of bispectra and trispectra from the Lifshitz scalar are generated because the action 
is not severely constrained by symmetries. As we discuss in our previous paper [41], however, only the local type 
bispectrum stems from the term which does not have the shift symmetry. In this paper, since we want to distinguish 
between the models where the equilateral type bispectrum are generated, we concentrate on the model with the shift 
symmetry. Moreover, only in the case with the dynamical critical exponent z — 3 we can have the scale invariant 
power spectrum. Thus we choose z = 3. 

The terms up to fourth order of the action of Lifshitz scalar with z — 3 and with the shift symmetry in the ultra 
violet is written as 



S=\j dtd 3 xa(t) 3 



+ j^*** + ^ {MA^HW + b 2 (^f) 



(CI) 



where M, bi, b 2 , t\, t 2 and £3 are constants. Then, the interaction Hamiltonian conforms with the non-linear terms 
of the action given by Eq. (CI). 

We can decompose the contributions of the trispectrum into six parts which are proportional to b\, 6162, b 2 , t\, t 2 
and £3. The scalar exchange contributions are the first three and the others are the contact interaction contributions. 
According to our previous paper [41], the shape functions of these contributions can be written as 

ah{se,ij) = M 4 (fcf + fcf + fcf + fcf + kf 2 )Mk 1 ,k 2 )f J (ks,k 4 ) 

r 2 3 ^k 1 k 2 k 3 k i k 2 2 (k 3 1 +kl +fc? 2 )(fcf + k\ + kl 2 ){kl + kl + kl + klY { ' 

q h(d,i) M 4 k 12 r l (k 1 ,k 2 ,k 3} k 4 ) 

r ~ 2 3 ^k 1 k 2 k 3 k i (k 3 + k 3 + k 3 + k 3 i )' 1 ' 

/i(fci, kj) — bi(k® + kj + kfj — kfk 4 — kfkj — kfkfj — kfkfj — k^kfj — k^k^), (C4) 

f 2 (k i ,k j ) = 6b 2 ^kpl, (C5) 

r^fci.fcj.fcfc.fe,) = 4t! {klk 2 3 {k k ■ ki) + k 2 k kf(ki ■ kj)) , (C6) 

r 2 (fc i ,fc J -,fc fc ,fc i ) = At 2 ((ki ■ kj) 2 (k k ■ k t ) + (k k ■ h) 2 (ki ■ kj)) , (C7) 

r 3 (ki,kj,k k ,ki) = 2t 3 ^(fcj • kj)(ki ■ k k )(ki ■ k t ) + (kj ■ ki)(kj ■ k k )(kj ■ k t ) 

+(k k ■ ki)(k k ■ kj)(k k ■ ki) + (ki ■ ki)(k t ■ kj)(k t ■ fe fc )) , (C8) 

where sty se ' % ^ and S^ cl '^ are shape function propotional to bibj and t{, respectively. Here, except for S^ cl ' 3 \ we 
can choose the reduced trispectrum such that it depends only on the 5 parameters k\ 1 k 2 , k%, k^ and ki 2 as is the 
case in DBI inflation (see Appendix A). Only S^ c% & cannot be reduced in such a manner because of the presence 

of the product of (fci • k 2 ) and (fci • k 3 ), etc. Therefore, in the case of S^ cl ' 3 \ we define the redused trispectrum so 
that it is symmetric under the transpose among k 1} k 2 , fc 3 and k 4 . 
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